Colloidal crystal growth at externally imposed nucleation clusters 
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We study the conditions under which and how an imposed cluster of fixed colloidal particles 
at prescribed positions triggers crystal nucleation from a metastable colloidal fluid. Dynamical 
density functional theory of freezing and Brownian dynamics simulations are applied to a two- 
dimensional colloidal system with dipolar interactions. The externally imposed nucleation clusters 
involve colloidal particles either on a rhombic lattice or along two linear arrays separated by a gap. 
Crystal growth occurs after the peaks of the nucleation cluster have first relaxed to a cutout of the 
stable bulk crystal. 

PACS numbers: 82.70.Dd, 81.10.Aj, 64.70.Dv, 81.10.-h 



While important steps towards a quantitative under- 
standing of homogeneous crystal nucleation out of the 
melt have been made in the past decade (for recent re- 
views, see [TJ [2]), work on the molecular principles of 
heterogeneous nucleation is still at its infancy [SI HI E] ■ 
Colloidal suspensions have served as excellent model sys- 
tems for nucleation, since the crystallization process is 
typically much slower than in their molecular counter- 
parts and the critical nucleus can be detected in real 
space [BJ. By using external fields, e.g., optical tweez- 
ers, it is possible to fix a cluster of colloidal particles 
and watch directly its impact on the rest of the colloidal 
suspension. If the crystal phase is slightly off-coexistence 
and the fluid is stable, it is possible to generate crystalline 
layers around such a cluster [71 03 IH] ■ 

In this Letter we study crystal growth processes into 
a metastable fluid. A cluster of fixed colloidal particles, 
which could act as a seed for heterogeneous crystal nucle- 
ation is arranged within the metastable melt. Whereas 
in homogeneous nucleation such clusters spontaneously 
form by means of thermal fluctuations, here they are ex- 
ternally imposed. We investigate whether they serve as 
initiators of crystal growth processes. Our motivation 
for this study is twofold: first, by imposing a seed cluster 
one can steer the crystallization behavior. Second, offer- 
ing complex cluster structures could lead to unexpected 
dynamical scenarios of crystal growth. 

We approach the problem using classical density func- 
tional theory (DFT) of freezing which is a microscopic 
approach to crystallization [T0J [TT]. DFT can be ex- 
tended to describe dynamics in strongly inhomogeneous 
Brownian fluids [12j [T3l [T4J [15] . Here it is put forward 
as the first full microscopic approach to the dynamics of 
crystallization. Our DFT results are backed by Brown- 
ian dynamics (BD) computer simulations. In principle, 
the dynamical DFT is superior to phase-field crystal the- 
ories of nucleation |16j . which operate on more coarse- 
grained length and time-scales and need phcnomcnologi- 
cal mobilities as an input. Therefore our results provide 
benchmark data to test the validity of more approximate 



theories. 

In detail, we study a model for a two-dimensional sus- 
pension of superparamagnetic colloids, exposed to an ex- 
ternal magnetic field which tunes their parallel dipole 
moments (TTJ [TB]. By using additional fields, such as 
optical tweezers, certain particles can be fixed in the sus- 
pension |T9l [20]. We first consider a stable fluid phase, 
realized for a weak magnetic field. In this fluid suspen- 
sion, colloidal particles are placed by optical tweezers 
into prescribed positions forming a cluster. Then the 
magnetic field is suddenly increased rendering the fluid 
metastable with respect to the stable hexagonal crystal 
and the tweezers are released. Two different kinds of clus- 
ter geometries are considered: In the first setup we study 
hexagonal clusters that are cut out of a perfect rhom- 
bic lattice while in the second setup two sets of linear 
crystalline arrays, separated by a gap, are examined. 

As a result, we observe that the kinetic pathway of 
the system is a two-stage dynamical process: first, on a 
sub-Brownian time scale, the peak positions of the exter- 
nally imposed nucleation cluster relax towards a cutout of 
the stable bulk crystal. Then, on a Brownian time scale, 
there are two further possibilities: either the relaxed clus- 
ter acts as a nucleation seed for further complete crystal 
growth or it dies out completely without stimulating fur- 
ther crystallization. Whether crystal growth occurs or 
not depends delicately on the compatibility of the initial 
cluster geometry with that of the stable bulk crystal in 
terms of strain energy. 

Our system is characterized by the pairwise inter- 
action potential u(r) = uo/r 3 , where uq is the inter- 
action strength. For the specific realization of two- 
dimensional paramagnetic colloids of susceptibility x ex ~ 
posed to a perpendicular magnetic field B, we have 
u o = (xB) 2 /2. The thermodynamics and structure 
depend only on one dimcnsionless coupling parameter 
r = uop 3 / 2 /hsT ', where p is the average one-particle den- 
sity and kgT is the thermal energy. 

It has been shown [TSJ, Q3J [14] that the static, classical 
DFT can be given an extension to dynamics to describe 
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overdamped, time-dependent, out-of-equilibrium systems 
in terms of a deterministic, time-dependent, and ensem- 
ble averaged one-particle density p(r, t). The time evolu- 
tion of p(r, t) is then governed by the continuity equation 



t<0: equilibrium liquid 
t>0: undercooled liquid 



dp(r,t) 
dt 



p(r,t)V 



SF[p(r,t)} 
Sp(r,t) 



(1) 



Here, D/ksT is the mobility coefficient originating from 
the solvent, ignoring hydrodynamic interactions. 

The equilibrium phase diagram of the system under 
study has been obtained using classical DFT [2T], which 
provides the intrinsic Hclmholtz free energy functional 
F[p(r)), a unique functional of the static one-particle 
density p(r) of the system. The functional F[p(r)] is 
minimized by the equilibrium one-particle density, where 
it takes the value of the system's intrinsic Hclmholtz 
free energy. The density functional is typically split 
into the ideal gas, an excess, and an external part, 
F [p(r)j = F id [p(r)j + F cx [p(r)] + F CKt [p(r)]. The ideal 
part is F id [p(r)\ = k B T J drp(r) {in [p(r)A 2 ] - l}, with 
A denoting the thermal de Broglie wavelength. F^ is of 
completely cntropic nature and leads to a simple diffusion 
term in Eq. (Q. The excess part F ex , originating from 
the correlations between the particles, is in this paper ap- 
proximated by the ansatz of Ramakrishnan and Yussouff 
to the DFT [22] . It is expanded up to second order in 
terms of density difference Ap — p(r) — p around a refer- 
ence fluid, where the fluid density p is chosen as the av- 
erage density of the inhomogeneous system: F ox [p(v)] ~ 
F ex (p) - \k B T fj drdr'Ap(r)Ap(r')4 2) (r-r';p). Here 

(2) 

Fcx(p) and c (r; p) are the excess free energy and the 
direct correlation function of the reference fluid of den- 
sity p, respectively. Finally, the external part is simply 
given by F ext = / drp(r)V(r), where V(r) is the 

total external potential. 

For both setups under study, the clusters of tagged 
particles are first, i.e., for times t < 0, held fixed in a 
thermodynamically stable, equilibrated fluid of density 
p at a coupling constant of T< = 10, which is well be- 
low the freezing transition at T ~ 35.7 [23 , obtained 
within the theory. For the equilibration of the fluid, 
Eq. (JT|) is numerically solved fixing the tagged particles 
by deep parabolic external potentials at the tagged parti- 
cle positions — in an experiment this could be achieved 
by using optical tweezers |19j . At time t = we turn 
the external pinning potential off and, at the same time, 
instantaneously quench the system to a coupling con- 
stant r> = 62.5, which is well above the freezing tran- 
sition and we observe the time evolution of the den- 
sity field for times (£/ t b) ^5 10> where tb = (pD) -1 
is the Brownian time scale. Eq. ([T| is numerically solved 
applying a finite difference method. The dimensions 
,(\/3/2)a of the rectangular peri- 
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odic box considered are chosen integer multiples n x , n y 
of the lattice spacing a = (2/y/3) 1 / 2 p- 1 / 2 of the perfectly 



FIG. 1: Sketch of the imposed, rhombic nucleation cluster of 
19 particles surrounded by a gray fluid. The angle 4> between 
the spanning basis vectors and the area of a unit cell A are 
also shown. 



ordered hexagonal crystal. 

The first setup under study comprises a rhombic nucle- 
ation seed of 19 tagged particles, arranged in a hexagon, 
as sketched in Fig.[l] The nucleus is characterized by the 
strain parameters A, the area of a unit cell which in the 
perfectly ordered hexagonal crystal equals A = l/p, and 
4>, the angle spanned by two of the nucleus axes. The 
size of the periodic rectangular box is 16a x 16(v / 3/2)a. 
In Fig. [2] snapshots of the time-evolving density field 
arc shown cxcmplarily for two clusters cut out from two 
compressed hexagonal crystals with parameters {Ap — 
0.7, cos 4> — 0.5) and (Ap — 0.6, cos <f) — 0.5), respectively, 
at times t/rs = 0,0.001,0.1, 1.0. While the former, less 
strongly compressed cluster grows into the equilibrium 
crystalline state, the latter collapses back into an under- 
cooled, metastable fluid within t/rs < 0.1. The growth 
dynamics of the stable nucleus is a two-stage process: In 
the first stage — on a sub-Brownian time scale t < 0.002 
- the positions of the seed's density peaks move to a 
cutout of the thermodynamically stable bulk crystal. In 
the second stage — on the Brownian time scale — the 
system crystallizes out of the relaxed cluster. 

Fig. [3] displays the "island" of growth in the (^4,cos<^)- 
parameter space, i.e., the set of parameters, for which 
the nucleus grows for t > 0. It is found that the "island" 
is nearly symmetric in cos 0, relative to the equilibrium 
value of cos (f> = 0.5 while it is asymmetric in unit cell area 
A about the ideal value of A = 1/p. This asymmetry is 
qualitatively validated by BD simulations [231 125] . 

Within the second setup we study the time evolution 
of a nucleation seed of two equal linear arrays along the 
y-direction, each comprising three infinite rows of hexag- 
onally crystalline particles, which are separated by a gap, 
as can be seen from the density map for t — in Fig. [4] 
These arrays, corresponding to an equilibrium crystal 
generated via a suitable external potential, are displaced 
relative to each other in y-direction by half a lattice spac- 
ing Ay — a/2. In between the two crystalline arrays there 
is an empty stripe of width Ax = V3a corresponding to 
one missing row of crystalline particles. In contrast to 
the first setup, the second setup corresponds to a config- 
uration with a huge local, non-affine strain relative to a 
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FIG. 2: Snapshots of the central region of the dimension- 
less density field p(r, t)/p of two colloidal clusters with strain 
parameters Ap — 0.7, cos <f> — 0.5 (left panel) and Ap = 
0.6, cos = 0.5 (right panel) at times t/r B = 0,0.001,0.1,1.0 
(from top to bottom; t/r B = 1.0 only for Ap = 0.7) [24] . 
Note that the images display only the system's central region 
of dimensions L x /2 x L y /2. 



perfect cut-out of a bulk crystal due to the gap. 

In order to keep the gap free of particles during the 
equilibration of the surrounding fluid for times t < 
we employ an additional strong external potential in the 
region of the gap. The dimensions of the periodic box 
within which Eq. (fTl) is solved numerically are now given 
by L x x L y = 64(y3/2)a x a. Snapshots of the central 
region of the density field p (r, t) are shown in Fig. [4] for 
times t/rs = 0,0.01,0.1,0.63,1.0 after the quench. 

Again, a two-stage dynamical scenario is observed: On 
a sub-Brownian time-scale of about 0.02 tb, the positions 
of the peaks drift to those of a perfect cut-out of the 
stable bulk crystal. This leads to a rapid filling of the gap. 
Then crystallization occurs on a Brownian time-scale. In 
Fig. [5] we plot the distances Xi(t) of the three crystalline 
density peaks and the distance of the crystal front xj(t) 
with respect to the center of the gap as a function of 
time. The latter is taken as the inflection point of the 
envelope function of the y-averaged density field. The 
theoretical curves are compared to BD simulation data 
of the same setup |26j obtained by averaging over the 
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FIG. 3: Stability "island" of the imposed nucleation cluster 
of 19 particles according to Fig. [T] The shaded region sep- 
arates the growth from the no-growth situation. The (blue) 
stars display the according boundaries for fixed cos<j> = 0.5 
obtained from BD computer simulation. The hexagon sym- 
bols indicate the way the seeds are deformed in the different 
regions of the parameter space. The (red) dots indicate the 
configurations for the snapshots in Fig. [2] 



y 







<: :3 




















































£88888888 



I -1 5 -1 -0.5 0-5 1 1.5 2 

• log, ip(rj/p) 

FIG. 4: Snapshots of the central region of the dimensionless 
density field p(r, t)/p of a linear nucleus of two times three 
infinite rows of hexagonally crystalline particles, separated 
by a gap, at times t/r B = 0,0.01,0.1,0.63,1.0 (from top to 
bottom) [24] . Note that the images display twice the system's 
central region of dimensions L x /4 x 2L y for better visibility. 



particle positions of 24.000 independent simulation runs. 
The two-stage picture is clearly confirmed. 

In conclusion, we have investigated by dynamical den- 
sity functional theory whether and how an externally 
imposed cluster of fixed particles acts as a nucleation 
seed for crystal growth if the particles are released and 
the system is quenched instantaneously from a stable 
to a metastable bulk fluid. If the imposed cluster is 
not too much strained relative to a cut-out of the sta- 
ble bulk crystal, it induces global crystallization. The 
kinetic pathway of the imposed cluster exhibits a two- 
stage scenario: the cluster structure first relaxes towards 
an appropriate cut-out of the bulk crystal before further 
growth. This two-stage process is unexpected since it is 
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FIG. 5: Time evolution of the distance of the linear array's 
three density peaks Xi(t) and of the crystal front Xf(t) with re- 
spect to the center of the gap as a function of time. Dynamical 
density functional theory results (lines) are compared against 
Brownian dynamics simulation data (symbols; the dashed line 
connecting the crosses is a guide to the eye). The arrows in- 
dicate the typical time scales on which the relaxation of the 
Xi is occurring and on which the crystal growth sets in, re- 
spectively. 



reversed in larger clusters which contain quite a large por- 
tion of the stable bulk crystal. In the latter case crystal 
growth starts at the edges but the inner elastic distortion 
anneals on a much larger time scale. For higher under- 
coolings, i.e., larger r>, the size of the stability island 
(Fig. 3) increases. 

Our predictions can be verified by real-space exper- 
iments on two-dimensional superparamagnetic colloidal 
particles confined to the air-water interface in an external 
magnetic field [TTltlH] ■ Qualitatively similar scenarios are 
expected for different repulsive interactions and in three 
spatial dimensions, which arc relevant for nuclcation and 
growth experiments in sterically and charge stabilized 
suspensions [HJ [37J HE! US] • In three dimensions, one may 
even induce the growth of metastable crystals and quasi- 
crystals imposed by suitable nucleation seeds [3D] . 
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